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Abstract 

We show that the ttNN coupling constant extracted model- inde- 
pendently from pp charge exchange is subject to a systematic correc- 
tion, and, more importantly, that the strong absorption in the critical 
region prevents a determination of the coupling constant to high pre- 
cision using this process. This attenuates the possible conflict with 
the value determined from the np charge exchange cross sections. 

The value of the vrNN coupling constant has become a topic of hot debate 
in the last few years In addition from the 'classical' determinations from 
ttN scattering and forward NN dispersion relations with a 'high' value of 
near 14.3, the Nijmegen group has argued a lower value of about 13.6 on the 
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basis of their extensive analysis of NN data using the Nijmegen model [0, 0]. 
The issue of its precise value is an important question not only because it is 
a fundamental constant to nuclear physics , but also because it is of great 
importance also for the understanding of chiral symmetry breaking. Thus, 
the experimental error in the pion-nucleon coupling constant is the princi- 
pal obstacle in the precise testing of the validity of the Goldberger-Treiman 
relation as predicted from chiral symmetry breaking This relation 

would give = 12.81, if it were exact, which is not expected, however. To 
resolve the problem unambiguously we have undertaken a detailed study of 
the possibility of using methods of extrapolation to the pion pole in NN scat- 
tering for an accurate and model independent determination of the coupling 
constant directly from data. A first report on this work, based on high pre- 
cision absolutely normalized data in the pole dominated region of np elastic 
charge exchange at 162 MeV, points to a high value of = 14.60 ± 0.30 [^]. 
The dominant source of uncertainty is the experimental overall normaliza- 
tion. New, as yet unpublished, additional data give a slightly lower and more 
accurate value g"^ = 14.4±0.2 10] !§[• The pole extrapolation method was orig- 
inally suggested by Chew using a polynomial expansion . Its practical use 
is fraught with problems associated with systematic errors and instabilities 
which must be controlled and well understood. We have tentatively solved 
this problem on the one hand by stabilizing the extrapolation procedure using 
appropriate comparison functions, on the other hand by exploring systemat- 
ical errors analyzing large numbers of 'pseudo experiments' generated from 
models with known coupling constants. This program is still being pursued 
on a wider body of np charge exchange data. 

The corresponding information is in principle also available from the pp 
nn charge exchange reaction, which has the same pole structure. Recently 
precise data on this reaction was obtained in the experiment PS206 by Birsa 



et al. [10| at 601 MeV/c (= 176 MeV kinetic energy) and the quality of 



these data are comparable and even superior to those of the best data on 



np charge exchange in a wide range of similar energies [11]. A preliminary 



analysis demonstrated strikingly that indeed both pp and np charge exchange 



data extrapolate approximately to the same pole residue ]11]. However, the 



further analysis [|T^] based on the classical Chew method gave the very low 
value /2 = 0.0708±0.0016±0.0011 (^^ = 12.80±0.29±0.20) (statistical and 
normalization errors, respectively) , about 10% below our value above and 
also well below the Nijmegen value. This value is a cause of concern and raises 
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the question of consistency with np charge exchange and even more seriously 
of the overall validity of the approach as such. The analysis did not attempt 
to analyze and evaluate the systematic errors in the procedure. The present 
note aims to elucidate these questions as well as to clarify the structure of 
the contributions to antiproton charge exchange and their relations to the 
pion pole terms. 

We first examine the contributions to the unpolarized differential cross 
section so as to illustrate the information that must be described by any 
extrapolation procedure. Here the total amplitudes and cross sections are 
defined in the usual fashion in terms of the five amplitudes a, b, c, d, e allowed 



by the invariance properties ||T3|, . Both the CM unpolarized cross section 
as well as the polarization transfer one are incoherent combinations of 5 
amplitudes: 

^(t) = l/2(|ap + |6p+|cp + |dr + |en; 

^(t) (i-ir„oon) = icr+Mr, (1) 

where t is the squared 4-momentum transfer from the proton to the neutron. 
Taking into account the difference between the proton (Mp) and neutron 
(Mn) masses one has in the centre-of- mass system for the proton momentum 
k and scattering angle 6: 



t = +2k^k^ - M2 + M2 cos e-2k' + M^- M;. (2) 

There are five regularized pion Born amplitudes and in this case the r- 
space 5-function has been subtracted [O] |T6|: 



«7r + &7r = 0; 

- 6^ = F(t); 
c^ + d^ = F{t){l-m{t))- 
c^-d^ = F{t); 

Cn = 0, (3) 



where 
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In Eq. (4) the charged pion mass is denoted by m,r and the range of the form 
factor is chosen to be A = 800 MeV. This corresponds to an rms radius of .6 
fm for the nucleonic pion source. 

It is in practice convenient to chose a representation for the amplitudes 
such that the pion pole term appears in one single term, namely [c + d)/ \/2 
for the pp case and {h-d)/^ioi the np one. This regroups the pole term 
completely into the real part of these amplitudes, respectively, and they 
combine to an exact square in the cross section. These are therefore crucial 
to the extrapolation and they are also the dominant terms in the physical 
region. All the remainder is represented by terms which slowly and smoothly 
decrease with increasing —t. The differential cross section and polarization 
transfer then become 

^(t) = l/4[|a + 6p + \a- ftp + |c + dp + |c - rfp + 2|ep)]; 

^(t) {l-K^mn) = \[\c + d\' + \c-d\'] (5) 

Let us now examine the behavior of the contributions to the cross section 
and their variation with t in more detail. Such a study is most readily done 
using a model which represents the main physics even if it only qualitatively 
reproduces the data. The Paris antiproton model is sufficient for this pur- 



pose and readily reproducible fT^- The characteristic shape and magnitude 
of the terms in eq. (5) are illustrated for the pion Born terms of eqs. (3-4) 
and for the Paris model in Figs, la and lb. Not unexpectedly, absorption is 
a prominent feature for antiproton charge exchange. Compared to np charge 
exchange the cross section is quenched by a factor 4-5 for momentum trans- 
fers of a few pion masses, the region most sensitive to the pion pole. This 
feature is strikingly apparent also in the comparison of the np and the pp 
pole extrapolation in fig. 2. There is an important difference in the physics 
between the two cases. In the np case the amplitudes have imaginary parts 
only from unitarity. They are mostly nearly real and there are in practice 
only two important non-pole terms. In the antiproton charge exchange real 
and imaginary non-pole amplitudes are both important and four such terms 
contribute substantially. Fortunately our method permits us to ignore the 
detailed dynamics of these terms. In both case the term containing the pole 
amplitude has a zero at t = — in the Born approximation independent 
of form factors. The corresponding minimum survives with a slight shift in 
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the model cross sections as seen in Figs, la and lb. Although the value 
is quenched in the antiproton case this suggests that similar extrapolation 
methods can, at least in principle, be used in both cases. Practice appears to 
confirm this. The minimum is accompanied by a steep rise of the differential 
cross section at small angles. It is this feature that is responsible for the 
overall minimum in the cross section, which is partly masked by slowly vary- 
ing background terms. In the case of np scattering this minimum transforms 
into a shoulder in the cross section, but the corresponding structure is also 
present in that case. 

We note in passing that polarization transfer in principle allow a clean 
separation of the terms containing the pion pole contributions from other 
terms as is apparent from eq. (5). However, the use of this possibility still 
requires the accurate knowledge of the unpolarized cross sections for a direct 
extrapolation to the pion pole. Since the eliminated terms are not dominant, 
even a perfect knowledge of the polarization is only a minor constraint on 
the pion coupling constant. The polarization information is of no particular 
use in the present context. 

We have first followed a procedure identical to that of Bradamante et 
al., analyzing their data using the Chew method |ll|], [12|. In this case one 



defines the function y(x), which extrapolates the experimental data smoothly 
to the pion pole: 

SX^ dcr '^"^ 

Here s is the square of the total energy and x = —t + rri^. We will use the 
charged pion mass 171^^ = 1 as the mass scale in the discussion. At the pion 
pole X = the Chew function gives 

2/(0) = ao = Ng'/gj,, (7) 

where N is the normalization of the experimental data, which may differ 
from the true value. The model independent determination of the coupling 
constant requires accurate single-energy data with absolute normalization 
N of the unpolarized differential cross section. The error in the coupling 
constant determined by such extrapolation methods is always proportional 
to TiV. 

We accurately reproduce the Bradamante analysis for the pp reaction 
using the Chew method and the results are given in Table 1. However, we 
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find minor discrepancies in their simultaneous analysis of the Uppsala 96 
MeV np scattering data, which they use as a comparison. We traced this to 
their arbitrary omission of three data points. Therefore, at this first step our 
analysis reproduces their low value of the coupling constant. 

The next problem is that of systematical errors in the extrapolation pro- 
cedure and true extrapolation errors. Here the pp case differs importantly 
from the np, since absorption is a major feature. Examination of the Chew 
function in fig. 2 shows that it has a pronounced minimum near 0°, the phys- 
ical point closest to the value to be determined. The extrapolation point is 
about a factor of 4 larger than this value and lies opposite to the main 
trend of the function to be extrapolated. As a matter of principle this is a 
particularly unfavorable situation in any accurate extrapolation procedure. 
Another way of stating this problem is to note that the absorption produces 
an important reduction of the pion pole effects in the physical region, which 
is to say that it reduces the sensitivity to the coupling constant which we 
want to determine. Even so there remains a substantial sensitivity on the 
pion-nucleon coupling constant as already demonstrated by Bradamante et 



al. [|Tl|, |l^]. The accuracy to which it can be determined is however a more 



delicate and quantitative question. 

To investigate this matter we first use models with known g"^ for which 
we can generate randomly 'data points' at the exact angles and with the 
same statistical errors as in the actual experiment. For this we generate 
10.000 equivalent pseudo experiments which we analyze in the same way as 
the actual experiment. As models we use a solution of a preliminary analysis 



of the Nijmegen group and the Paris potential prediction [|T3- However 
any other model with known pion exchange would do as well for the present 
purpose if it approximately reproduces the data. From Table 1 we find that 
the Chew extrapolation approach requires 6 polynomial terms to describe 
these 'experimental' data appropriately. From the results we conclude that 
the Chew extrapolation procedure systematically underpredicts the coupling 
constant by 2% {dg"^ ~ 0.3) for the case of 47 data points (n=6). This shift 
would appear even if the precision of the data were higher than at present. 
It is thus necessary to correct for such systematic shifts in the determination 
of the coupling constant using the Chew method. 

However, very little information on the pion pole is contained in the 
region of large momentum transfers. We therefore also truncated the data at 
Xmax = 4.9m^ (the 30 first points), so as to be able to use fewer parameters 
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in the extrapolation. One easily persuades oneself that the tensor amplitude 
makes it physically unreasonable to expect a description of the data with a 
polynomial of less than 5 terms in the original Chew polynomial extrapolation 
procedure and that data with large momentum transfer call for at least 6 
terms. 

In the case of 30 points, the analysis of the pseudodata for both the Paris 
and Nijmegen models in Table 2 demonstrates that 5 terms in the polynomial 
give a perfect fit to data such as these, although with a systematic downward 
shift of by 0.3 (Paris) to 0.5 (Nijmegen). The extrapolated value is very 
low (12.26) and it is not more than about 12.8 ± 0.5, even when corrected 
for the systematic shift. The shift disappears with one additional term in 
the polynomial, but at the cost of a much larger extrapolation error. For 
the corresponding experimental 30 data points the x^/DoF is only 0.46, 
which is an unexpectedly low value statistically. It reflects partly that the 
errors given by PS 206 consist of statistical and systematic errors added in 
quadrature ll^. Corrected for systematic shifts the Chew model in this case 
gives 12.84 ± .46 for n = b and 13.24 ±1.15 for n = 6. The low raises the 
possibility of some correlations in the systematic errors, in particular at low 
momentum transfers. In every case the systematic error of 1.5%, i.e. ±.2, due 
to the normalization uncertainty must be added to the overall uncertainty. 
These results are compatible with the Chew analysis in the larger data range 
with 47 data points. The latter appears to be more accurate with a value 

= 13.1±.3as seen in Table 1 for n = 6. This is somewhat disconcerting 
since the additional region is insensitive to the pion pole and the polynomial 
expansion is questionable over such a large region. 

We now turn to the question of the precision that can be attained reliably 
in the analysis, since the errors given above are only the formal statistical 
uncertainty of the Chew method. To investigate this point it is convenient 
to use the Ashmore method as described in ref. 0]. We expect this method 
to work, since both the background terms and the interference term vary 
similarly to the np case as already discussed above. It is thus the systematics 
of this method that will determine the precision which can be achieved. 
The Ashmore method parametrizes the differential cross section in terms of 
the regularized pion Born amplitudes, but has in addition a contribution 
simulating a p meson pole with adjustable strength and shape, described 
with a form factor and polynomial terms in x 0. This method permits 
a phenomenological simulation of the absorption effects, but the dynamical 
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description as such is not realistic. The virtue of the method is its ease of 
apphcation and that it permits the exphcit exploration of the influence of 
changes in the pion-nucleon coupling constant. We note from Table 2 that 
the Paris and Nijmegen model pseudodata are described in the Ashmore 
model for n=5 with / DoF=1.00 with a moderate systematic shift of 0.21 
and 0.29, respectively, for 30 data points. An additional parameter as used 
with n = 6 will thus overparametrize the analysis. 

In order to test the uncertainty of the coupling constant determined from 
the full range of data in the Ashmore model we fixed the coupling constant for 
n = 6 to 13.21, 14.10 and 14.43 and determined the corresponding values 
40.2, 51.0 and 58.6 for the 47 data points (41 degrees of freedom). This should 
be compared to the corresponding Nijmegen description for g"^ = 13.23, which 
is 48.6. Although the value of 13.23 gives a better in agreement with the 
other methods all of these values give a good representation of the data as 
is apparent from Fig. 3. We note in particular that the first few points near 
t = have a large contribution to the x^- They play thus an essential role in 
the discussion. For example, one half of the difference between the of the 
Nijmegen model with a coupling constant of 13.23 and the Ashmore model 
with 14.43 comes from the two points nearest t = 0. Since these points are 
also the ones with the largest systematic errors ^TO), they are an important 
potential source of systematic error in the extrapolation. If they are omitted, 
the value of increases by .1 to .2. Interestingly, all of the Ashmore curves 
have a generically similar shape which differs systematically from the one in 
the Nijmegen description on the level of a few %. These observations suggest 
that the value of g^ deduced from antiproton charge exchange cannot be 
considered to fully model independent at a precision higher than 0.5 to 0.7 
units (3 to 5%). 

In the case of np charge exchange, which we previously investigated, it was 
possible to improve the procedure considerably using a Difference Method 
[^]. This method relies on an extrapolation of the difference of the Chew 
function for data and a model with a known coupling constant. While this 
method is effective in the case of np charge exchange, it does not improve 
significantly on the previous methods in the present case. The reason is the 
scarcity of high precision information on the pp interaction other than charge 
exchange. This has the consequence that either the model relies heavily on 
the present data on charge exchange or lacks sufficient precision in the de- 
scription of higher momentum transfer to be useful. In the former case the 
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method would as expected accurately reproduce the of the model with 
considerable precision, but the argument would be circular. The method has 
the additional advantage of visually bringing out the details of the extrap- 
olation in the low t region at the level of precision under discussion here. 
It is therefore instructive to apply the method to the Paris model, which 
has qualitative agreement with data and incorporates constraints at large 
momentum transfer. 

The results of the extrapolation of the difference function yuix) — y{x) 
versus x are given in Table 2 for 30 data points and are displayed in figure 4. 
We first calibrated the method using pseudodata from the Nijmegen model. 
These demonstrate that such data can be perfectly described by n = 5 with 
X^/DoF = 1.00 and that the systematic extrapolation error is then negligible. 
However, applied to the actual data we find no improvement in the accuracy 
to which the coupling is determined. On the other hand Fig. 4 demonstrates 
that a major background has been removed. It is now easy to visualize 
the consequences of different values for the number of polynomial terms n 
in the extrapolation. A good description is obtained already for n = 4 
but with a substantial systematic correction. The values deduced for the 
coupling constant after the correction for the systematic shifts are in this 
case 13. 10 ±.16, 12. 65 ±.45 and 13. 24 ±1.15 for n = 4, 5 and 6, respectively. 
These extrapolations give similar results to the previous ones and they are 
mutually consistent. We have also applied the method to the full set of data. 
In this case (n = 6) is needed, but the conclusions do not change. 

We have already remarked that the absolute normalization of the cross 
section is a crucial number for the extraction of the coupling constant. Let us 
now discuss the uncertainties from this source, which is independent of the 
errors arising from the extrapolation uncertainty. The Nijmegen group has 
achieved a good description of the present data with /{data) = 1.035 with 
an integrated charge exchange cross section atot = 12.14 mb. The normal- 
ization error is most likely the experimental uncertainty in the experiment 
PS206 in view of the close fit to it. For the Paris model the corresponding 
cxtot = 13.45 mb. This means automatically that if this normalization were 
used to normalize the same experimental data, then the deduced coupling 
constant would increase by 5.3%, that is would increase by 0.7. This 
emphasizes the importance of this question. 

Some additional information on this point comes from other experiments. 
R.P. Hamilton et al. [jl9[ have performed a dedicated absolute measurement 
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of the integrated charge exchange cross section with energy in small steps. 
By interpolation between 596 and 608 MeV/c this corresponds to a value 
of atot = 11.80(11) at 601 MeV/c. The systematic error is stated to be less 
than 3 to 5 %, depending on the energy region. Using 4% as a reasonable 
estimate at the present energy this means a systematic uncertainty of about 
0.47 mb, completely compatible with both the experiment PS 206 [|1^ and 
with the Nijmegen description. Experimental information can also be found 
in Nakamura et al. pO[. In this case the emphasis was on the angular shape 



and the forward interference dip and not on normalization. Their integrated 
cross sections are larger than in the other experiments by 15 to 20%, but the 
systematic uncertainty was also larger (8.5%). In view of the circumstances 
we, as the authors, do not believe these larger values to be significant. In 
conclusion, there appears to be no obvious reasons to question the presently 
quoted normalization in the experiment PS 206, but we emphasize that this 
issue must be kept in mind as a potential source of problems. 

We have here critically examined the accuracy to which the ttNN cou- 
pling constant can be extracted from the recent precision data on pp charge 
exchange and the importance of systematic theoretical corrections to pole 
extrapolation procedures. Our conclusion is that the situation is less favor- 
able than for the corresponding np charge exchange reaction. This is due to 
the prominent role of absorption which reduces sensitivity to the pion in the 
most critical region of momentum transfers. In particular, it is not profitable 
to apply the Difference Method in the present context. Here, the accuracy 
does not increase and it was this method that was the key to high precision 
for np charge exchange. In spite of this limitation and the more important 
role of systematics in the extrapolation procedure it is possible to extract 
the coupling constant to a good degree of precision, though with an errors 
of about 4 to 5%. Depending on the detailed procedure we find a range of 
plausible deduced values from 12.8 to 13.2 from the data. The statistical 
errors are typically of order .45 using only the range of data sensitive to the 
pion information. Formally, a higher statistical precision is achieved using 
a larger range of data. However, we found in model studies that a nearly 
indistinguishable description of the data is achieved with values for the cou- 
pling as large as 14.1. We therefore believe that caution should be used in 
quoting the formal errors above. We recommend a value = 13.0 ± .7 with 
an additional systematic error of 0.3 from the overall cross section normal- 
ization. The value that is extracted in this way is a direct determination 
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the coupling constant for the ttNN system. This value is low, but in view 
of the uncertainties it is compatible with the value for the coupling constant 
deduced from other sources. There is no obvious discrepancy with the value 
deduced from np charge exchange. 

We are indebted to F. Bradamante and A. Martin for details of the anal- 
ysis of the antiproton data PS 206, M. Lacombc for the solutions to the Paris 
antiproton potential, J. Blomgren and N. Olsson for discussions concerning 
on np charge exchange data and to R. Timmermans for information on the 
results of the Nijmegen analysis of the antiproton charge exchange data PS 
206. 
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Table 1: Results applying the Chew method to the Birsa et al. data points 
at 176 McV [10] and to the corresponding Nijmcgcn and Paris antinucleon 
pseudo-data (see text). The model coupling constants are gjsiijmegen — 13.23 
and glparis — 14.43 with bg^ the systematic shift from the true model value. 



Chew Method 




PS206 


'Nijmegen' 


'Paris' 


n 


x'iDoF 


xVDoF 


g' 


5g' 


xVDoF 


g' 


5g' 


5 


1.28 11.78 ±0.15 


1.77 


11.53 ±0.16 


1.70 


2.54 


12.42 ±0.14 


2.01 


6 


0.90 12.76 ±0.27 


1.00 


12.90 ±0.27 


0.33 


1.01 


14.17 ±0.24 


0.26 


7 


0.79 11.67 ±0.58 


1.00 


13.10 ±0.51 


0.13 


1.00 


14.47 ±0.45 


-0.03 
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Table 2: Results applying different extrapolation methods to the 30 first 
data points of PS 206 [10] as well as to the corresponding Nijmegen and 
Paris model pseudodata. 



Chew Method 




PS206 


'Nijmegen' 


'Paris' 


n 


xVDoF 




xVDoF 




Sg' 


xVDoF 


g' 


Sg' 


5 


0.47 


12.26 ±0.46 


1.00 


12.85 ± 0.44 


0.48 


1.00 


14.11 ±0.39 


0.32 


6 


0.46 


13.12 ± 1.15 


1.00 


13.11 ± 1.17 


0.12 


1.00 


14.38 ± 1.03 


0.05 


Ashmore Method 


n 


r/DoF 


9 


X'/DoF 


9 




X'/DoF 


y 




4 


0.44 


12.40 ±0.23 


1.04 


12.70 ±0.22 


0.53 


1.08 


13.90 ±0.20 


0.53 


5 


0.46 


12.43 ±0.34 


1.00 


12.94 ±0.33 


0.29 


1.00 


14.22 ±0.30 


0.21 


6 


0.46 


13.43 ± 1.55 


1.00 


13.14 ± 1.61 


0.09 


1.00 


14.38 ± 1.43 


0.04 



Difference Method 




Paris-PS206 


Paris-'Nijmegen' 


n 




oF (/ 




oF 


9 




3 


4.71 


15.13 ±0.06 


3.67 




14.98 ±0.06 


-1.75 


4 


0.70 


13.69 ±0.16 


1.10 




13.82 ±0.16 


-0.59 


5 


0.45 


12.62 ± 0.45 


1.00 




13.20 ±0.43 


0.03 


6 


0.46 


13.03 ± 1.15 


1.00 




13.02 ± 1.17 


0.21 
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Paris model contributions to pp ^ nn at 176 MeV 
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Figure 1: The magnitude and structure of the different contributions in eq. 
5 to the differential pp fin charge exchange cross section at 176 MeV: a) 
using the pion Born amplitudes of eqs. 3 and 4; b) using the Paris antinucleon 

J_1 rri1__ 1_ -1- -1 I _ i Jl2 




^ ' ' ' ' ^ ' ' ' ' L 

2 4 

X = (-t+l) (pion mass units) 



Figure 2: Comparison of the Chew pole extrapolation for np and pp charge 
exchange at 162 McV [4] and 176 MeV [10], respectively. The reference 
couphng constant is 14.1. For details, see text. 
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Figure 3: Fit to the data of Birsa et al. [10] by the Nijmegen group 
{gljq^ = 13.23) [18] and using for the n = 6 Ashmore model with 3 different 
values for the nNN coupling {g"^ = 13.23, 14.10 and 14.43, included in the 
definition of n = 6). Note the similar quality of the descriptions. 
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Figure 4: Extrapolation to the pole of the difference between the Chew- 
function for the Paris model [17] and for PS206 [10] at 176 MeV for n=3, 4 
and 5. Note that n = 4 or 5 is required by the data. 
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